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ABSTRACT 



We investigate the properties of the most optically faint sources in the GOODS-N area (R AB > 26.5). Such extremely optically faint 
populations present an uncharted territory despite the fact that they represent an appreciable fraction of the X-ray sources in the 
GOODS-N field. The optically faint sources are believed to contain either red AGN at moderate redshifts or possibly QSO at very 
high redshift. We compile our sample by first finding the 3.6/jm IRAC counterparts of the X-ray sources and in turn by searching 
for the optical counterparts of the IRAC sources. 35 sources do not have counterparts in the W-band Subaru optical images. Of these, 
18 have HST ACS counterparts while the remaining have no optical counterparts. The vast majority of our 35 sources are classified 
as Extremely Red Objects (EROs) on the basis of their V^oe - K s lower limits. Their derived photometric redshifts show that these 
populate moderate redshifts (median z ~ 2.8), being at markedly different redshifts from the already spectroscopically identified 
population which peaks at z ~ 0.7. The Spitzer IRAC mid-IR colours of the sources which have no HST counterparts tend to lie 
within the mid-IR colour diagram AGN "wedge", suggesting either QSO, ULIRG (Mrk231) templates or early-type galaxy templates 
at z > 3. A large fraction of our sources (17/35), regardless of whether they have HST counterparts, can be classified as mid-IR 
bright/optically faint sources (Dust Obscured Galaxies) a class of sources which is believed to include many heavily absorbed AGN. 
The co-added X-ray spectrum of the optically faint sources is very flat having a spectral index of T as 0.87, significantly flatter than 
the spectrum of the X-ray background. The optically faint (R > 26.5) X-ray sources constitute more than 50 per cent of the total X-ray 
population at redshifts z > 2 bearing important implications for the luminosity function and its evolution; considering X-ray sources 
with 2 < z < 4 we find good agreement with a modified Pure Luminosity Evolution (PLE) model. 
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1. Introduction 

X-ray surveys provide the most efficient method for detecting 
AGN (Brandt & Hasinger 2005). This is because X-ray wave- 
lengths can penetrate large amount of interstellar gas and re- 
veal the AGN even in very obscured systems. The deepest X- 
ray surveys to date detect a large number of sources down to 
a flux of ~ 2 x lQ-^ergcm-V 1 in t he (0.5-2.0) keV band 
dAlexander et al.L l2Q03h ILuo et all 120081) . The vast majority of 
them are AGN dBauer et all 2004 ) with a surface density of 
about 5000 sources per square degree. Optical foll ow-up obser- 
vation s have identified a l arge fraction of them dBarger et al., 
I2003at ICapak et all 12004 iTrouille et all [2008), revealing that 
the peak of the redshift distribution based on spectroscopic iden- 
tifications is at z — 0.7. However, a large number of the X-ray 
sources remain optically unidentified hampering our understand- 
ing of their nature. In particular, a large fraction of faint X- 
ray sources (~ 50%) lacks a spectroscopic identification (e.g. 
ILuo et all l2010h and the redshi ft estimate is made with pho- 
tometric techniques. Still. lAird et al.l {2010) estimate that about 
one third of the X-ray sources in the CDFs do not have optical 
counterparts down to Rab ~ 26.5. 

The nature of these optically faint sources remains puz- 
zling. Two interesting scenaria have been proposed to explain 
their nature. First, the very faint optical emission could be 
the result of copious dust absorption; in the Chandra deep 
fields the majority of high f x /f sources shows clear evidence 



of obscuration in their individua l and stacked X-ray spectra 
dCivano. Comastri & Brusal 120051). a result which confirms pre- 
vious findings ( Alexander et alll200lh . On the other hand, a very 
faint, or the lack of a n optical counterpart cou ld indicate a very 
high redshift source ( Koeke moer et all 120041) . In this case, the 
reason of a high X-ray to optical ratio is that the optical bands are 
probing bluer rest-frame wavelengths, which are more obscured 
or intrinsically fainter if they fall blueward from the Lyman 
break. At the same time the observed X-ray wavelengths cor- 
respond to high ene rgy rest-frame wavel engths which are less 
prone to absorption. iLehmer et alj d2005l) used a Lyman break 
technique to select high redshift galaxies in the Chandra deep 
fields and found 1 1 B435, V(,o6, and (775 dropouts among the X- 
ray sources, with possible redshifts z > 4. 

A crucial diagnostic for the nature of optically faint galax- 
ies is their infrared emission. The optical and ultra-violet light 
which is absorbed is re-emitted in infrared wavelengths, there- 
fore a high infrared to optical r atio can be us e d as a criterion of 
high obscuration. For example, iHouck et al.l d2005l) discovered 
a number of sour ces with extremel y high 24ytzm to optical lumi- 
nosities (see also lDaddi et alll2007h . These sources, nicknamed 
DOGs (Dust Obsc ured Galaxies) are located at 1.5 < z < 2 
dPope et all 120081) . A large fraction of them is probably as- 
sociated with Compton-fhick QSOs dFioreetall I2OO8I [20091: 
iTreister et al., 2009; Georgantopo ulos et alll2008l) . 

In this paper we explore the properties of optically faint 
(^ab <: 26.5) X-ray sources in the GOODS-N area. Previous 
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studies of some of the optic ally faint source s in th is field have 
been performed in the past; Alex ander et al.l (1200 l l ) studied the 
"brigh test" optically faint sources (with / > 24) and lCivano et all 
(2005) investigated the properties of Extreme X-ray Objects 
(EXOs), having an X-ray to optical flux ratio > 10. Instead, 
our sample focuses not only on EXOs but on all optically faint 
sources. Additionally, we take advantage of the full 2 Ms X- 
ray exposures, but more importantly we make use of the mid- 
IR Spitzer observations in order to find the true counterparts 
of the X-ray sources. We perform anew the search for optical 
identifications of the X-ray sources in two steps. First, we cross- 
correlate the X-ray positions with the Spitzer IRAC 3.6 pm cata- 
logue using a likelihood ratio technique. Then we cross-correlate 
the IRAC positions with the ICapak et al.l d2004l) 7? -band optical 
observations. We define the optically faint sources as those with 
Rab ^ 26.5, i.e. beyond the flux limit of this catalogue. Then 
we look for optical counterparts in the fainter HST observations. 
After deriving photometric redshifts (section 4) we study their 
X-ray (section 5) as well as their mid-IR properties (section 6). 



2. X-ray, optical and IR data and counterpart 
association 



We use the 2 Ms CDFN catalogue of lAlexander et al.l (§003), 

- 1 in the 0.5-2.0 keV 

,-2„-l 



with a sensitivity of 2.5 x 10 17 erg cm 2 s 



band and 1.4 x lO^ergcirr^ 1 in the 2.0-8.0keV band. The 
infrared (Spi tzer) data come from the Spitzer-GOODS legacy 
programme (Dicki nson! 120041) . This includes observations in 
the mid-infrared with IRAC and MIPS which cover most of 
the CDFN area. Typical sensitivities of these observations 
are 0.3 pJy and 80/Jy for the IRAC-3.6yum and MIPS-24//m 
bands respectively. Optical coverage of the GOODS-North 
field is available both with ground-based (Subaru-SuprimeCam) 
and space observations (HS T-ACS). The Suba ru sensitivity is 
26.5 mag(AB) in the R band (ICaoaket all 12004]) and t he HST is 
27.8 mag(AB) in the zsso band dGiavalisco et al., 2004). We also 
make use of near-infrared (Ks -band) observations with Subaru- 
MOIR CS, reaching a limit of 23.8mag(AB) dBundv et all 
l2009h . 

The X-ray catalogue (lAlexander et all 120031) has 503 
sources detected in one of the hard (2.0-10.0keV) soft (0.5- 
2.0keV), or full (0.5-10.0keV) bands. Of these sources, 
348 fall into the region covered by the IRAC observations. 
We make an initial simple cross-correlation of the positions 
of the sources of the two catalogues to check their rela- 
tive astrometry. We find good agreement in the RA axis 
(< RA(Xray) - RA(IR) >= -0.05, cr = 0.31), but there is 
a significant difference in the astrometry in the DEC axis 
(< DEC(Xray) - DEC(IR) >= -0.29, <x = 0.35, see Fig.Q); we 
update the positions of the IR sources accordingly before pro- 
ceeding to the search of IR counterpart s to the X-ray sources. 

W e use "likelihood ratio" method (Suth erland & Saunders! 
1992) to find the counterparts: we calculate the likelihood ratio 
(LR) of an infrared source being a real counterpart of an X-ray 
source as: 



LR = 



q(m)f(x,y) 
n(m) 



where q(m) is the expected magnitude distribution of real coun- 
terparts, f(x,y) is the probability distribution function of posi- 
tional errors, a gaussian in this case, and n(m) is the magnitude 
distribution of background objects. We calculate the q(m) func- 
tion by subtracting the magnitude distribution of the background 
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Fig. 1. Distance in RA and DEC between X-ray and infrared 
sources. There is a significant shift in declination between the 
two catalogues which we correct for before making the correla- 
tion. 
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Fig. 2. Same as Fig.[TJbut for the Spitzer (without including the 
correction introduced to match the X-ray astrometry) and optical 
catalogues. The shift in declination is again corrected for before 
correlating. 

(IR) sources we expect to find within the search radius near the 
target (X-ray) sources from the magnitude distribution of all the 
initial counterparts. The normalization of q(m) is done using 



/Mb,, 
oo 



q(m)dm = Q(M hm ) 



where Q(M\ lm ) is the probability that the IR counterpart is 
brighter than the magnitude limit Mi; m ); in practice it is the final 
fraction of X-ray sources with an IR counterpart. Given that the 
positional offsets of the X-ray sources in the CDFN at larg e off- 
axis angles can be as large as 1"— 2" ( lAlexander et alll2003l) . we 
use a large initial search radius for the X-ray - IR counterparts 
(4"). 

The choice of an optimum likelihood ratio cutoff is a com- 
promise between the reliability of the final sample and its com- 
pleteness. The reliability of a possible counterpart is defined as: 



Ri = 



LRj 



£LRj + (l-Q(M llm )) 



where j refers to the different IR counterparts to a specific X-ray 
source. We chose the likelihood ratio cutoff which maximizes 
the sum of completeness and reliab i lity o f the matched catalogue 
in the same manner as iLuo et all (l2010h . The completeness is 
defined as the ratio of the sum of the reliabilities of the counter- 
parts with LR > LR\ im with the number of the X-ray sources in 
the area mapped by IRAC, and the reliability is defined as the 
mean reliability of counterparts with LR > LR\i m . In the X-ray - 
IR case LRii m = 0.15, which gives a reliability of 99.2%. 

We detect 330 IR counterparts with LR > 0.15, we then op- 
tically inspect the positions on the IRAC images of the X-ray 
sources lacking a counterpart and find 12 cases where the IR 
source is clearly visible but blended with a nearby source in a 
way that the source extraction algorithm could not distinguish 
them. Adding these cases, the final number of counterparts is 
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342 out of the 348 X-ray sources in the IRAC area. Assuming 
that 0.8% of the counterparts are spurious, the final efficiency is 
97.5%. 

In order to look for optical counterparts to the X-ray sources 
we use the positions of their IR counterparts. The positional ac- 
curacy of IRAC is better than that of Chandra, moreover the NIR 
is more efficient in detecting AGN than the optical emission, due 
to the reprocessing of the ionizing radiation of the AGN through 
circumnuclear and interstellar dust to IR wavelengths. As a re- 
sult, the efficiency of IRAC in finding counterparts for the X-ray 
sources ( 97.5%) is higher than this of the optical survey (78.1% 
if we follow the procedure described above for the Subaru-/? cat- 
alogue), in spite of the optical being deeper than IRAC (see also 
iLuoetall HoH for the CDFS case). 

The optical catalogue of ICapak et all d2004l) has 47450 
sources detected in the R band, and of those 14763 fall into 
the area sampled by IRAC. Before searching for optical counter- 
parts, we check again the relative astrometry of the optical and 
(uncorrected) IR catalogues. We find again a similar astrometry 
difference in the DEC axis (< DEC(opt) - DEC(IR) >= -0.22, 
cr = 0.20, Fig|2]), while the RA positions are well within the er- 
ror (< RA(opt) - RA(IR) >= 0.05, <x = 0.17). We again correct 
the IR catalogue to match the optical positions. 

We use the likelihood ratio method to select the optical coun- 
terparts to the IR sources. Given the smaller PSF of the optical 
images we use a smaller initial search radius (3") and follow 
the same procedure to select the optimum likelihood ratio limit. 
With LR\ lm - 0.25 which this time gives a mean reliability of 
98.4% we find optical counterparts for 8437 of the 10595 IRAC 
sources. The recovery rate is 78.0%. We perform also an inde- 
pendent search for optical counterparts of X-ray sources which 
do not have an infrared counterpart, or the latter is a merged 
source as described above, and find a secure optical counter- 
part for 12 such cases; 1 1/12 of the "confused" IR sources have 
an optical source related to the X-ray position, while 1/6 of the 
IRAC non-detections is detected in optical. Finally, we optically 
inspect the Subaru R-band images of X-ray source lacking an 
optical counterpart and find that in 17 the optical counterpart 
is clearly vi s ible b ut with no identification in the catalogue of 
ICapak et~ai1 d2004t) . Two common reasons for that are that the 
source is saturated, and therefore it does not have reliable pho- 
tometry, and that it lies close to a very bright source and therefore 
it is not detected by the source extracting algorithm. 

The X-ray, optical, and infrared properties of all the 
sources of the common CDFN-IRAC area can be found in 
ftp .mpe .mpg . de/people/erovilos/CDFNJR.OPT/ 

3. Sample selection 

In this study we are interested in X-ray sources which are too 
faint t o be detected by typ ical ground-based optical surveys, like 
that of Ca pak et"aT]d2004|) . Such sources represent a sizable frac- 
tion of the overall X-ray population (~ 15%). They are gen- 
erally not covered with spectroscopic surveys, which are typi- 
cally magnitude limited with Rn m < 26.5 and often do not even 
have photometric redshi fts. Such cases have been studied be - 
fore in detail (e.g. lAlexander et all 120011 : iMainieri et al.L 12005): 
here we approach them using their infrared properties and X-ray 
spectra. 

There are 310 X-ray sources within the common area of 
all the surveys used in this study (2 Ms CDFN - Subaru- 
Suprime(optical) - Subaru-MOIRCS(///T) - GOODS-ACS - 
GOODS-IRAC - GOODS-MIPS). Of these 310 sources, 42 are 
too faint to be detected in the /?-band by Subaru-Suprime and are 



assumed to have R^b > 26.5. No /^-detected CDFN source has 
Rab > 26.5. Because in this study we are based on the infrared 
properties to examine the nature of optically faint sources, we 
exclude 7 of the 42 sources from our studied sample for the fol- 
lowing reasons: 5 are not detected with IRAC, and 2 are blended 
with nearby sources, so their infrared photometry is not reliable. 
Our final sample consists of 35 sources, whose X-ray, optical, 
and infrared properties are listed in Tab.Q] Among these sources 
we expect a mean number of 0.35 spurious encounters (practi- 
cally none), as the mean reliability of the IRAC counterparts for 
these 35 cases is 99.0%. 

For these 35 sources we search the K$ catalogue of 
iBundv et alJ ((2009) for counterparts. These Subaru-MOIRCS 
images cover an area slightly smaller than IRAC and we limit 
our field to the common area. As the number of sources is 
small, we look for Ks -band counterparts by eye; we can this 
way easily distinguish any Kg -band sources being physically re- 
lated to a nearby IRAC source. We find a K$ counterpart for 
29/35 sources; their distances from the IRAC sources are all 
< 1 .4", while for the 6 non-detections the distance of the near- 
est Ks source is always > 3.9". We t hen search the HK '-band 
images from the UH-2.2m telescope (Capa k et all 120 04) to de- 
tect any sources bright enough in the HK' band but not detected 
in R. Using sextractor (Berti n & Arnoutsl [T996). we identify as 
a source 4 adjacent pixels with fluxes above 1.2 times the lo- 
cal background rms. We find an HK' detection for 19 of the 
sources in our sample. The Subaru-MOIRCS area has very deep 
optical imaging wit h the HST-ACS as part of the GOODS sur- 
vey. The catalogues dGiavalisco et all [20041) are based on detec- 
tions in the zsso band and are publicly available. We search the 
sources of our sample for HST counterparts using the likelihood 
ratio method, as described in the previous section with an initial 
search radius of 1 .5 arcsec. We use this method because the PSFs 
of Spitzer-IRAC and HST-ACS are different within a large factor 
and there could be multiple ACS sources within one IRAC beam. 
18/35 sources of our sample have an HST-ACS counterpart. 

We note here that som e of the sources in Ta b.Q] appear with 
full optical photometry in iBarger et alJ d2003al) . some of them 
even with a photometric redshift. In cases of X-ray sources lack- 
ing an optical counterpart, like the ones in Tab.Q] Barger et al. 
d2003al) measured the optical fluxes directly from the Subaru im- 
ages using a 3"diameter aperture centered on the position of the 
X-ray source. In doing so there is a high probability that light 
from a neighboring source enters the aperture, moreover center- 
ing on the X-ray positions causes a loss in positional accuracy, 
which is essential when measuring the flux of a "non-visible" 
source. Fig. [3] shows the R -band images of the 35 sources of 
Tab.Q] with 4"-radii circles on the X-ray positions (as large as 
the initial search radius) as well as contours representing the 
IRAC flux. We can see that in some cases the X-ray and IRAC 
positions differ significantly and that the 1.5"aperture centered 
on the X-ray position would often be contaminated by nearby 
optical sources. 

The sample presented in Tab.Q] has both similarities 
and differences wit h the high X-ray to optical ratio sample 
dCivano et all 120051). and the optically faint sample presented 
by [Alexander et al. (2001). In Fig. [4] we plot the R magnitudes 
of the optical counterparts against the (0.5-8.0) keV flux as filled 
circles and the sources in our sample (all lower limits) as crosses. 
The shaded area is the lo g(./x/./opt) = 0_± 1 area where the bulk of 
AGN are expected (e.g. Elvi s et allll994l) . Our sources have, on 
average, a higher / x // op t value than the bulk of the AGN popula- 
tion, but co uld not be securely consi dered as high / x // op t (with a 
value > 1 0. iKoekemoer et all 120041) . The HST V(,ob magnitudes 
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provide a clearer picture. In Fig. [5] we plot the magnitudes 
against the X-ray flux keeping the same symbols as in Fig.|4] 
The V606 magnitudes for non-detected sources are assumed to 
be V606 > 27.8, the same as the detection threshold in the Z5&0 
band. The true limits however are likely to be higher, as the typ- 
ical V606 - Z 850 colour of AGN h osts is > 0.5 and increasing 
with redshift (Sanc hez et al.L l2004): the 2<x limits of sources 246 
and 317 which are detected in the zsso band and not in the ^606 
band are 29.6 and 29.9 respectively. In Fig. [5] we see that the 
sources of Tab.Q] are clearly relatively hi gh f x /f np t, all havin g 
l°g(/x/ /opt) > 0.52. However, compared to lCivano et al.l d2005). 
this is neither fully a high / x // op t sample, nor a complete one, as 
many sources with log( f x / f n pt ) > 1 are not inclu ded. Moreover, 
the high / x // pt sample of ICivano et al.l d2005l) includes much 
brighter sources, as bright as R ~ 23. 

Compared to lAlexander et alj d2001l). our s ample probes 
much fainter sources: the lAlexander et al. (2001) sample has a 



cut at / = 24. As a result, their sources have relatively low X-ray 
to optical flux ratios, in many cases log(/ x // opt ) < 0. 

4. Photometric redshifts 

We u se the EAZY code (Bramme r. van Dokkum & C oppi. 
120081) to calculate photometric redshifts. We used 4 ACS (B435, 
^606; his, and zsso) optical, HK' (when available), K$ and 2 
IRAC bands (3.6 fjm, 4.5 ymi) to constrain the photometric red- 
shifts. We do not use the 5.8 fim and 8.0 /ww bands because they 
are more sensitive to the properties of (interstellar and circum- 
nuclear in AGN cases like here) dust, something which would 
add extra parameter s which would have to be con sidered in the 
template fitting (see Rowan-Robin son et all [2008). 

The results are also shown in Tab.[T] the HST magnitudes 
shown as lower limit are detections with a < 2cr confidence in the 
respective bands, and the 2<x flux is used to calculate the mag- 
nitude. The photometric redshifts based on these lower limits 
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Fig. 4. Optical versus X-ray flux for all the X-ray sources in the 
common area (see text). Crosses mark the sources of Tab. Q] and 
the shaded area marks log(/ x // opt ) = 0+1 where the bulk of 
AGN are expected. 
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Fig. 5. Same as Fig.gjbut for the HST- ACS V m optical flux. 
The V606 > 27.8 upper limit for non detection is likely to be an 
under-estimation (see text). 



are considered less reliable. Note that source 369, despite being 
very faint in the optical, has a repo rted spectroscopic redshift of 
z — 2.914 in lChapman et alj 12005). The position of the spectro- 
sc opic source com es from the radio (VLA - 1.4 GHz) catalogue 
of lRichardsl (l2000h and agrees within 0.15"with the IRAC posi- 
tion. The spectrum is typical of an AGN (CIV line) and the X-ray 
position is 1.1 "away, so we are confident that this is the correct 
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redshift 

Fig. 6. Redshift distributions of CDFN sources in the common 
CDFN-GOODS-MOIRCS area. The black histogram represents 
all optically detected sources with the shaded histogram repre- 
senting spectroscopic redshifts. The red histogram (refering to 
the right axis scale) represents the optically faint sources. 



counterpart. The photometric redshift derived by EAZY (2.80) 
is close to the spectroscopic value. 

For sources with no ACS detection, it is very challenging 
to derive a photometric redshift based only on the IRAC and Ks 
bands. In these cases we are forced to use all IRAC bands, which 
adds the extra variable of the dust properties. A good example of 
poor photo-z accuracy in cases where only m id-infrared fluxes 
are fitted can be found in lSalvato et al.l d2009l) . The redshift de- 
rived is typically 2 < z < 3, but the constrain is weak and the 
broad-band spectrum in some cases can be equally well fit with 
SEDs shifted to z > 3.2. The reason for that is that in cases where 
the IRAC SEDs are monotonic in f v the fit cannot be easily con- 
strained. A good photo-z estimate with a value of z ~ 2.5 is 
derived in cases of a blue [5.8]-[8.0] coloufl which can be fitted 
with the red-NIR bump of a galactic ( no AGN) SEP due t o mod- 
erate temperature interstellar gas (see lSalvato et al.Ll2009l Figure 
8). From Tab.[Q we can see that of the 14 sources with full IRAC 
photometry and no ACS detection, 5 have blue [5.8]-[8.0] colour 
(98, 140, 156, 198, and 372). The median photometric redshift 
for these five sources as derived by EAZY is z — 2.39, whereas 
the median photo-z for the red [5.8]-[8.0] colours is z — 4.40. 
The two samples define different populations in terms of pho- 
tometric redshifts with a confidence level of 99.5 %, according 
to a K-S test. Moreover, the '7/5T-detections" and "HST-non- 
detections" populations are also different within 98.9 %, having 
median redshifts 2.61 and 3.48 respectively. 

The redshift distribution of the sources of Tab.Q~]is compared 
with the distribution of all X-ray sources in the common CDFN- 
GOODS-MOIRCS in Fig.|6] The black histogram represents all 
optically detected sources and the shaded histogram represents 
so urces with sp ectroscopic redshi fts, taken from the catalogues 
of lBarger etakl d2003al) and iTrouille et all d2008l) . The red his- 
togram (refering to the right axis scale) represents the sources of 
Tab.Q] The two distributions are remarkedly different, with the 
sources of Tab.Q]being at higher redshifts. 



1 We define here blue as f v5 8/zm > f v go pm , corresponding to [5.8] - 
[8.0] < 0.64 
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5. X-ray spectral analysis 

The data have been analyzed using the CIAO v.4. 1 analysis soft- 
ware. The source spectra are extracted from circular regions with 
variable radii so as to include at least 90 per cent of the source 
photons in all off-axis angles. There are 20 X-ray observations 
comprising the total ~2 Ms exposure. We extract the source 
spectrum and auxiliary files for each observation separately us- 
ing the CIAO SPECEXTRACT script. Then we use MATHPHA 
ADDARF and ADDRMF tasks of FTOOLS to merge the spec- 
tral products for each source. The background files are calcu- 
lated from source free regions in each observation and are again 
merged using the MATHPHA tool. 

We explore the X-ray properties of the 35 sources in our sam- 
ple using the XSPEC vl2.5 package to perform X-ray spectral 
fittings. For the sources with adequate count statistics (net source 
counts > 200, 5 sources), we use the x 2 statistic technique. The 
data are grouped to give a minimum of 1 5 counts per bin to en- 
sure that Gaussian statistics apply. We adopt an absorbed power- 
law model and attempt to constrain the intrinsic absorption col- 
umn density A^h (i-e-, having subtracted the Galactic absorption) 
and the power-law photon index F. 

For the sources with limited photon s tatist i cs (ne t counts 
<200), we use the C-statistic technique (ICashl 1 19791) specifi- 
cally developed to extract spectral information from data with 
a low signal-to-noise ratio. In this case, the data are grouped to 
give a minimum of 1 count per bin to avoid zero count bins. 
We try to constrain the intrinsic column densities using an ab- 
sorbed power-law model with F fixed to 1.8. In both cases, the 
spectral fittings are performed in the 0.3-8 keV energy band 
where the sensitivity of the Chandra detector is the highest. 
The estimated errors correspond to the 90 per cent confidence 
level. In Tab. [2] we present the spectral fitting results for the 35 
sources comprising our final dataset. Source 107 has a very flat 
photon index, which is likely to be the result of a r eflection- 
dominated spectrum (see Georgantopoulos et al., 2009, for a de- 
tailed analysis of the candidate Compton-thick sources in the 
CDFN). In this case the hydrogen column density in Tab.|2]is not 
correct, the true A^h of the source being > 10 24 cirT 2 . There are 
two more Compton-thick sources (199 and 369) in our sample 
which are transmission dominated. In these the X-ray absorp- 
tion turnover is redshifted at low energies (5-6 keV) and an 
observer's frame A^h ~ 3 x 10 23 cirT 2 is measured, which yields 
a rest-frame A^h > 10 24 ctrT 2 when the (1 +z) 265 correction is ap- 
plied. The most secure Compton-thick source is 369 for which 
there is spectroscopic redshift available and thus a more reliable 
determination of the rest-frame column density. Note that this 
source is not included in the CDFN Compton-thick sample of 
iGeorgantopoulos et al.l (2009), because it has a flux lower than 
their adopted flux limit of 10~ 15 ergcirT 2 s _1 in the 2-10 keV 
band. 

Next, we compare the mean X-ray spectral indices by 
coadding the individual X-ray spectra in the total 0.3-8 ke V band 
for four different sets of sources (see Table [3). We fit a sim- 
ple power-law model to the data. It is evident that our sources 
present a very hard X-ray spectrum. The total spectrum of all 
35 sources has T « 0.9 which is harder than the coadded spec- 
trum of all t he det ected X-ray sources in this band (F = 1.4, 
iTozzi et alTl 12001 1) . Note that the deri ved mean spectru m is 
comparable to the spectrum d erived by iFiore et all d2008l) and 
iGeor gantonoul oset al.l (120081) for Dust-Obscured-Galaxies i.e. 
sources defined as having fiA^mlfR > 1000; the coadded spec- 
trum of DOGs in our sample has T = 1.05. The sources with 
IRAC-only detections are significantly harder than those with 



Table 3. Co-added X-ray spectral properties 



Sample 


No 




r 


Total 


35 


0.87 


±0.05 


HST 


18 


1.02 


±0.07 


IRAC-only 


17 


0.56 


±0.08 


DOGs 


17 


1.05 


±0.06 



optical HST counterparts. This difference cannot be attributed to 
the fact that they have on average different redshifts and hence 
different K-corrections, as it would work in opposite directions. 
The HST detected sources are intrinsically softer than the ones 
completely lacking an optical counterpart. 

6. Optical Properties 

In Fig|7] we plot the Vm - hi5 vs (775 - Z&50 colours of the 
18 sources in our sample which have an ACS detection. Red 
points mark the objects of our sample and black dots are HST 
sources associated with an X-ray source. All the HST detec- 
tions are shown in comparison in blue points. The open circle 
with the right arrow represents source 317, which has a low sig- 
nificance (< 2cr) detection in the and (775 bands and its 
Vgo6 - '775 colour cannot be accurately determined. The coloured 
lines in Fig. [7] trace the colours of various templates with red- 
shift with points marking z - 0, 1,2,3. The templates are the 
IColeman. Wu & Weedmanl ([l980) galaxy templates for ellipti- 
cal, spiral, and irregular galaxies. Dotted lines correspond to 
z > 4, where the blue edge of the V(,06 filter is redshifted to the 
Lyman break (912 A). The templates used are zeroed at shorter 
wavelengths than the Lyman break, so the dotted lines are ap- 
proximations of the colours at high redshift. The sources of our 
sample appear significantly redder than the overall X-ray popu- 
lation; 55.6% (10/18) of them have 2775 -Z&50 > 0.8 compared to 
14.4% (42/291) of the overall X-ray population; the probability 
that the /77s - z&sa colours of the red point being a random sam- 
ple of the black points in Fig[7] is < 0.01%. The colours of the 
red objects appear to be consistent with the elliptical template at 
Z = 2 ± 1, which is backed up by the photometric redshifts of 
Tab.Q] Red optical colours are oft en associated with early-type 
morphologies (see Bell et aUl2004l) . 

In Fig. [8] we plot the HST-ACS cutouts of the galaxies of 
Tab.Q] Each thumbnail is a combination of the four ACS (B— V— 
i-z) images to increase the signal-to-noise ratio and the contours 
represent the IRAC 3.6 ^m flux, as in Fig. [3] The crosses mark 
the positions of the HST sources with an inner radius of 0.5"and 
an outer radius of 2". We can see that for sources recovered with 
the HST the morphologies cannot be determined. 

Based on their near-infrared (Ks) to optical colours, many 
of the sources in our sample are as sociated with Extremely Red 
Objects (EROs; lElston et all 119881) . E ROs are usually defined 
as galaxies having R - K > 5 (see also lAlexander et al.L [2002, 
for a selection based on the /-band with / - K > 4). If we use 
the fl-band limit (R AB > 26.5 => /? Veg a > 26.3), 13/29 of the 
sources detected in Ks are EROs. However, the true optical flux 
of the sources is fainter than this limit; if we use the V606 ACS 
magnitude for HST detected sources and Vab > 27.8 => Vyega > 
27.7fl 28/29 of the ^ s detections have V 60 6 - K s > 5. Source 
335 has V 60 6 - K s > 4.86 and it is not detected with the HST, 

1 The detection limit in the zsso ACS band is 27.8(AB) and according 
to the colours of HST detections the V 606 limit is expected to be even 
higher. 
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Table 2. X-ray properties 



IAU A03 rest-frame Nn T Counts log(L2.o-io.okev) 







xl0 22 cm 2 






ergs 1 


CXO J123605.83+620838.0 


73 


< 1.56 


i -^-0.4 

1.8 


175 


44.00 


CXO J123608.61+621435.1 


84 


50 


20 


43.61 


CXO J123613.02+621224.1 


98 


AO 9 + 17.1 


1.8 


127 


44.15 


CXO J123614. 14+621017.7 


100 


121 3 +2L5 
15 +20 Z *1 


1.8 


56 


44.07 


CXO J123615. 83+621515.5 


107 


56 +0 41 

U ' J -0.31 

1.8 


211 


44.08 


CXO J123621.07+621303.3 


125 


8.5** 


33 


42.53 


CXO J123621.94+621603.8 


129 


1.8 


64 


44.03 


CXO J123623.66+621008.7 


140 


1.8 


69 


43.33 


CXO J123627.53+621218.0 


151 


26.8!** 


1.8 


52 


43.92 


CXO J123628.78+621 140.0 


156 


62.8!^ 5 


1.8 


42 


43.27 


CXO J123631.25+620957.3 


167 


< 16.0 


1.8 


35 


43.11 


CXO J123634.48+620941.8 


181 




1.8 


66 


43.45 


CXO J123636.90+621320.0 


196 


< 298 


1.8 


31 


43.86 


CXO J123637.26+620637.5 


198 


17.3!" 


1.8 


46 


42.91 


CXO J123637.28+621014.2 


199 


480!^ 


1.8 


26 


43.65 


CXO J123638.94+621041.5 


204 


13-2!l| 


1.8 


78 


43.46 


CXO J123639.65+620936.4 


208 


< 25.8 


1.8 


27 


42.88 


CXO J123642.11+621331.6 


220 


< 2.93 


1.8 


44 


43.01 


CXO J123647.94+621019.9 


246 


qq 7+12.4 


1.8 


118 


43.93 


CXO J123648.28+621456.2 


250 


1.60"$ *2 


9 O+0.7 

1.8 


170 


- 


CXO J123656.56+621513.1 


290 


138!ii 4 


81 


43.93 


CXO J123657.91+622128.6 


299 


7 05 +235 

* u -2 94 


1 7 1 +0.22 
i. ; i_ 23 

1.8 


787 


44.72 


CXO J123658.74+621459.2 


302 


< 159 


34 


42.97 


CXO J123659.32+621833.0 


307 


103.8! 27 8 { 


1.8 


67 


44.21 


CXO J123701.62+621 146.2 


317 


< 10.2 


1.8 


17 


42.99 


CXO J123702.43+621926.1 


321 


16.37! 2 ^ 

1 cn+1.50 #9 


i qo+0.22 

A ' 7Z, -0.28 

1 8 
1 .o 


433 


44.10 


L,A.W J 1ZJ /Uj.oJ + OZl J Ju.y 


DLy 


_o 




CXO J123705. 12+621634.8 


335 


67 3 +24 ° 


1.8 


125 


44.35 


CXO J123712.09+621211.3 


369 


1380+ 1 ' 40 


1.8 


31 


43.84 


CXO J123713.65+621545.2 


372 


1.8 


40 


44.38 


CXO J123713. 84+621826.2 


374 


31 9 + " 1 


2 16 +0 - 48 
1.8 


200 


43.03 


CXO J123725.50+621707.3 


417 


<6.3 


155 


43.58 


CXO J123733.98+621624.0 


434 


90!« 


1.8 


27 


43.44 


CXO J123737.04+621834.4 


445 


17.13«|] 


9 19+0.31 
^ -0.40 


492 


44.06 


CXO J123750.22+621359.3 


470 




1.8 


105 


44.16 



*1: Column density is more likely > 10 24 cm 2 (see text). 
*2: Observer's frame column density 



so it is highly probable that it is too an ERO. Of the six sources 
not detected in K 5 , 100, 151, 250, 329 and 374 are not detected 
with the HST either, so they could be associated with EROs and 
181 has V606 - < 5.03. Morphologically, EROs are mix of 
early and late-type systems dCimatti et all 12003b iGilbank et all 
120031) and the higher r edshift ones tend to be more late-type 
( Moustaka s et all |2004 . 

7. Infrared properties 

We plot the mid-infrared colours of the sources of Ta b.Q] in 
Figlll This diagram has been used by IStern et alj d2005) to se- 
lect AGN by their mid-infrared colours; a high fraction (up to 
90%) of broad-line AGN is located inside the "wedge" marked 
with the black solid line in Fig|9] The ffiiT-detected and non- 
detected sources of Tab.Q] are plotted with filled and open cir- 
cles ^j^snectrvel^^ The_coloured lines represent the colours of 
the Coleman et al. (1980) galaxy templates as well as a QSO 
SEP (lElvis et all Il994h for < z < 8, and points mark 
z = 0,1,2,3,4,5,6,7. 

The sources of Tab. [TJ are located both inside and outside the 
wedge in Fig|9] in the area where the red/optically-faint mid- 



IR AGN of iGeorgantopoulos etai] d2008l) lie. There is a sepa- 
ration in the colours of HST detections and non-detections, the 
latter population having mid-infrared colours mostly inside the 
wedge (10/14), and the rest being equally distributed (9/18). 
This difference reflects on their [5.8]-[8.0] colours, a K-S test 
shows that HST detections have bluer [5.8]-[8.0] colours than 
non-detections with 99.0% (2.6 cr) significance, and it is what 
causes the photometric redshif ts of HST non-detecti ons to be 
higher. The galaxy templates of Cole man et al.1 (Il980h enter the 
wedge at z — 3 so the photometric redshift of those sources is 
z > 3 if they are fitted with those templates. However, the QSO 
template is in the wedge independent of the redshift because of 
its power-law shape, in fact this is the reason why this diagram 
is a selection method of mid-infrared AGN. So, a red [5.8]-[8.0] 
colour and the position inside the wedge can be explained either 
with a QSO SED or a high redshift (z > 3) galaxy SED. 
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Fig. 8. HST cutouts of all sources listed in Tab.Q] The images are combinations of the four ACS bands available to increase the 
signal-to-noise ratio, the contours represent the IRAC flux, and the crosses are centred on the GOODS source positions (if any) with 
an inner radius of 0.5"and an outer radius of 2" 



8. Discussion 

8.1. High or intermediate-z sources ? 

In an earlier study of optically faint sources in the CDFN, 
lAlexander et al.l d2001l) concluded that they are moderately ob- 
scured AGN at redshifts z = 1-3, leaving a small margin for very 
high redshift QSOs in the optically fainter subsample. However, 
their modest optical magnitude cutoff (7 > 24) includes sources 
that have "normal" X-ray to optical ratios (/ x // op t ~ 1). Studying 
more extreme cases based on fainter op tical flux es or higher X- 
ray to optical ratios, iBarger et al.l (l2003bh and iKoekemoer et"ai1 
(2004) suggest that the existence of very high redshift objects 



cannot be definitely ruled out. In this study we compose a sam- 
ple of the optically faintest AGN (R > 26.5) with robust in- 
frared identifications which is 83.3% complete, including 35 of 
the 42 R > 26.5 X-ray sources in the common CDFN-GOODS- 
MOIRCS area used. 

8.1.1. Multi-/t investigation 

For the 18 sources with an HST identification we have calcu- 
lated their photometric redshifts using up to 8 optical and in- 
frared bands and the results show that they are indeed at moder- 
ate redshifts, with a median z = 2.61, while all have z < 3.65. 
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Fig. 7. Optical colours of HST detected sources. The optically 
faint sources of Tab.Q]appear in red and the X-ray sources appear 
in blac k. For comparison, all HST sources of iGiavalisco et all 
(2004) ar e plotted in blue. Colour lines track the colours of the 
IColeman etaT] dl980h SED templates with redshift. 
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Fig. 9. Infrared colours of the sources listed in Tab.Q] HST 
detected sources are plotted in filled circles and HST non- 
detections in open circles. All 3.6 //m detected sources are plot- 
ter in blue dots. Th e colour lines track the colours of the 
IColeman et all d!980t) SED templates with redshift and the black 
lines mark the reg ion where infrared-selected AGN are located 
(ISternetal.Ll2005h . 



The remaining 17 sources with no optical HST detection present 
larger difficulty in constraining their redshift, and we have at- 
tempted to derive photometric redshifts using the infrared bands. 
The redshifts for sources with mid-infrared colours outside the 
"wedge" of Fig.|9l or alternatively with f v5 . &flm > fv&.o^m (6 
sources) can be better constrained and they have a median red- 
shift of z = 2.39. We assume that in these cases the host galaxy 




1000 



24 tun 



(MJy) 



Fig. 10. 3.6 /im vs. 24/iin emission. Dots represent the overall 
3 .6 /urn population with a 24 fim counterpart, crosses the counter- 
parts of X-ray sources, black filled circles the sources of Tab.Q] 
detected with the HST and red open circles sources of Tab.Q] 
not detected with the HST. The lines mark log(/3 e^mlfiA^m) - 
0,-1,-1.5. 



dominates the mid-infrared colours which are well fitted with 
non-AGN SEDs at 2 < z < 3 (green, cyan, yellow, and magenta 
lines in Fig . |9j> . The 9 sources with no optical HST detection and 
fv5.%iim < /v8.0/im have mid-infrared colours compatible both 
with normal galaxy (or host galaxy) templates at high redshifts 
or with QSO templates without a strong redshift constrain (see 
§@); we try to constrain their redshifts using photometry from 
lower energy infrared bands (24 jjm, see lSoifer et all 12008) 

The position of the QSO templates inside the wedge is the re- 
sult of the heating of the circumnuclear dust with the radiation of 
the AGN, which results in a red power-law SED (f v oc y" with 
a < -0.5). Th e circumnuclear dust can also e xtingu ish the op- 
tical emission; iDunlop. Ciras uolo & McLurJ (120071) claim that 
optically faint sources can be well fit with dusty and extremely 
reddened (Ay — 4) SEDs at moderate redshifts (z ~ 2.5). We 
therefore check the mid-infrared 24 //m emission of the sources 
of our sample to examine the dust properties. Of the 17 sources 
with no HST detection, 7 are detected by MIPS, including 5/9 
sources with /vs.s^m < /vs.o/zm- In Fig.[l0]we plot the 24 /vm 
emission with respect to the 3.6 /im emission. As dots we plot 
the overall population (non-AGN), as crosses we plot the coun- 
terparts of the X-ray sources (AGN) and as black filled and red 
open circles we plot the optically faint sources with and without 
an HST detection respectively. We can see that the optically faint 
sources are in general fainter in 3.6 yum emission with respect to 
their 24 //m emission, which indicates dust extinction of optical 
and near-infrared wavelengths, unless there is a strong 24//m 
component from the AGN. The 24 fim emission upper limits of 
the HST non detections can in most cases be explained with dust 
obscuration, in cases where fc.e^m < 5 /vJy. The harder X-ray 
spectra of the HST non-detections also suggest that (despite their 
somewhat higher redshifts) these sources are more obscured, and 
this is what causes the faint optical fluxes. 
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The four sources with no HST detection, /s.s^m < /ko/an, and 
no 24 /mi detection (84, 151, 204, and 335) are possibly associ- 
ated with high redshift objects. However, the 24 /mi upper limits 
of three of them (151, 204, and 335) do not rule out dust obscu- 
ration given their low 3.6 /mi fluxes (< 3.6 /dy). Source 84 with 
/3.6/jm = 6.546 /dy has an fi.b^nl fiA^m ratio lower limit simi- 
lar to the overall AGN population with no clear signs of dust 
absorption and is therefore a high redshift candidate. Its pho- 
tometric redshift based on the 4 IRAC bands (4.50) is weakly 
constrained and not reliable. 



8.1.2. Dropouts 



The Lyman-break technique (e.g. ISteidel et all I2003I) is often 
used to select high-redshift objects, and is based on the sud- 
den drop in the flux of the broad-band spectrum at wavelengths 
smaller than 912 A. iLehmer et all (120051) used the HST obser- 
vations to search for dropout sources among the Z8 5 o-detected 
sources in the CDFN. They found two Ve,o6-dropouts at z > 5 
and one /775-dr opout at z > 6. One o f the two V606-dropouts 
(source 247 in Alexander et al., 2003) is not included in Tab.fT] 
because it is near a brighter optical source (1.6 arcsec separation) 
and the two are blended in the IRAC ima ge. It has a spectro - 
scopically confirmed redshift of z — 5.186 (Barg er et alll20 03a). 
The other V606-dropout (source 246) has a red (775 - Z&50 colour 
(1.1), like most of the sources in the ffiiT-detected sample, and 
its V606 - '775 colour is a lower limit (V6O6 - '775 > 1-275 if we 
take the V(,o(, 2 cr limit), being compatible both with early-type 
templates at z ~ 1 .5 and with very high redshift late-type tem- 
plates, depending on the true V(,o(, - z'775 colour. This source is 
also observed spectroscopical ly, but neither a redshift nor a spec- 
tral ty pe could be derived (Barger et al., 2003a; iTrouille et al.L 
2008). It has no MIPS 24 fxm detection, but its mid-infrared 
colours ([5.8] - [8.0] = 0.265) argue against a high redshift, they 
are explained by normal galaxy templates with 2 < z < 3 agree- 
ing with the photometric redshift (2.80) derived using 7 optical 
and infrared bands. 



The (775 -dropout of ILehmer et al.l (120051) is source 317. Its 
photometric redshift (3.25) is based in lower significance B435, 
V606, and (775 measurements and it is not reliable. This the source 
with the highest 24 /mi flux in our sample (/^m = 664 /dy), 
which would yield a mid-infrared luminosity in the order of 
10 26 W Hz 1 in the 3.5 /mi rest- frame band if the source was ly- 
ing at z ~ 6. Its low /3.6//m//24//m ratio however (10 -L1 ) implies 
dust obscuration according to Fig. [10] making it more likely to 
be optically faint as a result of obscurarion rather than high red- 
shift. 



The dropout technique is proposed to select high redshift ob- 
jects, it is based however in 2-3 bands, usually including upper 
limits. All sources in our sample are "optical dropouts" in the 
sense that they are not detected in wavelengths shorter than a 
limit, in this case the /?-band optical, and taking the HST obser- 
vations into account, almost half of them are z-dropouts. When 
considering their multi-band properties however we see that a 
small fraction of them (if any) can be high redshift {z > 6) 
sources, as there are other processes which can cause the faint 
optical fluxes. Especially since we are examining AGN-hosting 
systems the SEDs are a combination of the AGN and the host 
galaxy and a simple colour selection can be misleading. 



8.2. Obscured AGN ? 

The colour-colour diagram of Fig. [7] shows that the sources 
which are detected by the HST have on average redder optical 
colours than the overall X- ray population and are c ompatible 
with elliptical galaxy SEDs. | Rovilos & Geo rgantopo ulosld2007l) 
and iGeorgakakis et al.l ([2008 ) have shown t hat red-cloud AGN 
are obscured post-starburst systems. Instead, ferusa et al .1 (120091) 
attribute the red colours to dust reddening rather than an evolved 
stellar population. The fact however remains that X-ray AGN 
with red optical colours have a high fraction of obscured sources. 

A numb er of s tudies (e.g. Alexander etafl 1200 lb 
iMignoli et all l2004t iKoekemoer et all 120041) relate high 
X-ray-to-optical ratios to the Extremely Red Objects sample, 
with R - K > 5. The sources of Tab.[TJhave by definition high 
X-ray-to-optical ratios. Using the HST magnitudes, we get 
log(/x// op t) > 0.52 between the ACS-V band (A ef{ = 6060 A) 
and the (0.5-10)keV X-ray band, assuming an optical limit 
V606 > 27.8 (AB) for non detected sources (see §[3] and Pig.[5j. 
They also have very red colours with V606 - Ks > 5. X-ray 
detected EROs are assumed to be low-Eddington obscured 
AGN and they show on average high / x // op t values dBrusa et all 
2005), like our sources. The Z/5T-detected subsample which 
shows red optical c olours have s imilar characteristics with the 
high X/O EROs of lMignoli et al.l d2004l) . although this is a much 
fainter sample. Their sample is detected in the /^-band and 
has bulge-domin ated morphologies, d ominated by their host 
galaxies (see also Maiolin o et all 12006). 

The sources in our sample are selected to have faint op- 
tical fluxes, and thus naturally many of the m are associated 
with dust obscur ed galaxies (DOGs). DOGs dDev et all [2008; 
iPope et all [20081) are a class of objects, selected by their high 
mid-infrared to optical fluxes ( f v (24 um) / f v (R) > 10 00); they 
are found at a redshift of z — 2 (e.g. IPope et all 120081) and are 
thought to be marking a pha se of bulge and black- hole growth. 
Hydrodynamic simulations dNaravanan et alll2009l) suggest that 
DOGs are rapidly evolving from starburst to AGN-dominated 
systems through mergers in cases where / v (24/mi) > 300 /dy. 
Alternatively, in cases where / v (24/mi) < 300 /dy the evolu- 
tion is secular, that is through smaller gravitational pert urbations. 
Their morphological characteristics dMelbourne et alll2009l) re- 
veal less concentrated systems for lower luminosities. 

Considering the Veo6 magnitudes, and a Va)6 > 27.8 in cases 
of HST non-detections, all objects in Tab.fT] with a MIPS detec- 
tion (17/35 sources) have DOG characteristics^, and based on 
24 /mi upper limits the rest are not incompatible with DOGs, as 
the f v (24 fim) / f v (V(,06) upper limit is always > 1250. Therefore 
half of our sample have characteristics consistent with DOGs 
and this is only a lower limit. Only one source (317) has 
/ v (24//m) > 300 /dy. Their red optical colours (HST detec- 
tions) are suggestive of early-type bulge-do minated morpholo- 
gies dBell et all 12004; Mign oii et all |2004|) . although the mor- 
phologies cannot be determined from direct observations, in line 
with their low 24 /mi fluxes (see also lBussmann et aTLl2009l) . In 
conjunction with the obscured nature of the AGN, this means 
that the host galaxy is dominating the optical light. 

On the basis of X- r ay sta c king analysis in de ep X-ray 
fields, iFiore et al.l (I2008L 120091) : iTreister et al.1 (120091) propose 
that DOGs may be hosting a large fraction o f Compton-thick 
sources (see also Georgantopoulos et all 120081) . The sources of 
our sample appear to be more obscured than the overall X-ray 



3 They all have f v (24fim)/f v (V 606 ) > 1000, except source 196 with 
/v(24/tfn)// v (V 60 6) > 890 
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population, the a verage spectral ind ex being F = 0.87 com- 
pared to F = 1.4 dTozzi et al. L 1200 it) . The stacked spectrum of 
X-ray detected DOGs in the Chandra deep fields is F ^ 0.7 
dGeor gantonou los et all 120081) . comparable to these of our total 
and DOG samples. The average T of our DOG sample is very 
close to the value derived by Georgantopoulos et al. (2010, sub- 
mitted) using all DOGs in the CDFN regardless of their optical 
detection (or lack of it), a sample broadly overlapping with ours. 

According to Tab.[2l the average column density is a few 
times 10 23 cm~ 2 . A few Compton-thick AGN may be also present 
in our sa mple, this includes the refle ction dominated source 
(107; see Georgantopo ulos et al.L 2009). Two more "mildly" 
Compton-thick sources i.e. with column densities A^h ~ 10 24 
ctrT 2 are directly identified on the basis of their absorption 
turnover entering the Chandra passband. Finally, we note that 
a number of our sources, (10/35) of Tab.|2] have unobscured 
hard X-ray luminosities L x > 10 44 ergs~' and intrinsic A^h > 
10 23 ctrT 2 making them members of the QS02 class. It is inter- 
esting that these sources do not reveal their QSO nature except 
in X-rays. 

8.3. X-ray Luminosity Function Incompleteness 

In the CDFs a number of X-ray sources lack an optical coun- 
terpart, even in the deepest ground-based optical surveys. More 
specifically, within the common CDFN - GOODS - MOIRCS 
area there are 35 sources which are not detected at a magnitude 
limit of R ~ 26.5 and 17 sources not detected at Z&50 ~ 27.8. 
These populations represent 1 1.3% and 5.5% of the X-ray popu- 
lation detected in this area. These numbers should be considered 
as lower limits, as there are 7 more X-ray sources with no opti- 
cal counterpart not included in Tab.Qjbecause either they are not 
detected with Spitzer or have unreliable photometry. Moreover, 
there is a num ber of stars, normal galaxies, or ultraluminous 
X-ray sources ( Ho rnschemeier et all 120031: iBauer et all |2004) 
among the 310 X-ray sources, so the fraction of non optically 
identified X-ray AGN is even higher. 

This sample of optically faint X-ray sources is not random in 
their redshift distribution; the redshifts calculated in §|4]and the 
K-S test performed show that the optically faint X-ray sources 
have significantly higher redshifts than the overall population 
and the HST non detections even higher. The redshifts listed in 
Tab.Q]are in 34/35 cases higher than 1.5. If we count the sources 
within the comm on GOODS - IRAC - MOIRCS area with an op- 
tical detection in lCapak et alJ d2004l) and z > 1 .5 their number is 
41. So, the incompleteness of the X-ray z > 1.5 population with 
an optical identification and a redshift estimation is ~ 50% and 
it becomes even higher at higher redshifts; 29/35 of our sources 
have z > 2. 

This has implications in the calculation of the X-ray lu- 
minosity function at high redshift (z > 2 — 4) and its evolu- 
tion (e.g.lSilverman et aD.l2008HAird et alll2008l:lYencho et all 
2009; Air d et all 12010). More specifically, this incompleteness 
affects the faint end of the XLF, as its "knee" at z=2 is at 
logL x ~ 44.5 ergs" 1 (logL, = 44.4 - 45.0 ergs -1 depend- 
ing on the lumin osity evolution model, ISilverman et all 12008; 
I Aird et all 1201 Oh and 34/35 of the sources in Tab.|2]have unob- 
scured X-ray luminosities L x < 10 44 ' 5 ergs~' in the 2 - 8keV 
band. 

In order to test how the inclusion of the sources in Tab.Q]can 
affect the luminosity function at high redshift, we calculate the 
luminosity function of AGN in two redshift bins, 2 < z < 3, and 
3 < z < 4, w ith and without their c ontribution, using the method 
described bv lPage & Carreral (2000). We use the combination of 
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Fig. 11. The space density of AGN in the CDFN in two redshift 
bins (2 < z < 3, left, and 3 < z < 4, right) and two luminosity 
bins (10 43 < L2.o-8.okeV < 10 44 , and 10 44 < L2.o-8.OkeV < 10 45 ), 
with no correction applied for optical non-detections. Open sym- 
bols are based on the sources in the C DFN - GOOD S - MO IRCS 
area with a redshift determination in Barg er et"a"D (l2003al) and 
iTrouille et alJ d2008l) and filled symbols are based on the com- 
bination of those sources with those of Tab.Q] The solid and 
dashed lines rep resent the best-fi t mod- PLE and LDDE models 
respectively (see ISilverman et ail, 12008). 



the CD FN catalogues of iBarger etail d2003al) and ITrouille et all 
d2008l) and select sources in the combined CDFN - GOODS - 
MOIRCS area. We do not make any distinction between spec- 
troscopic and photometric redshifts and we do not correct for 
optically unidentified sources. This has an effect on the derived 
luminosity function only when we do not include the sources of 
Tab.Q] as the final sample is nearly complete. 

The result is plotted in Fig[TT] We can see that the inclu- 
sion of the new sources severely affects the calculated luminosity 
function. We note that in most cases in the literature the incom- 
pleteness caused by optically non detected sources is treated by, 
either taking an optically bright sample which is spectroscopi- 
cally complete, with the side effect of undersampling t he high 
redshift - low luminosity cases fe.g. lYencho et all 12009). or us- 
ing simplistic assumptions, like that their red shift distribution 
is dir ectly connected to their X-ray flux (e.g. ISilverman et all 
2008), which is generally not true, or treating the incomplete- 
ness as a free parameter when fitting the d atapoints to deter - 
mine the luminosity function evolution (e.g. lAird et all 12010). 
In FigQT| we also compare our dat apoints with models o f 
AGN evolution, which we adopt from Silverm an et alJ (|2008). 
These are a modified Pure Luminosity Evolution (mod-PLE) 
model (see Honkin Z Richards & Hernquistl 120071) (solid line) 
and a Luminosity Dependent Density Evolution (LDDE) model 
(dashed line). Our datapoints are remarkably close to the predic- 
tion of the mod-PLE model while the non-corrected points are 
closer to those of the LDDE model. A thorough investigation 
of the luminosity function evolution however would require tak- 
ing into account sources at all redshifts and luminosities which 
would vastly increase the reliability of the statistics. This anal- 
ysis is beyond the scope of this work, but we caution the fact 
that the optically faint sources, despite being a small minority of 
the X-ray sources can impact the LF in high redshifts because of 
their redshift distribution. 
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9. Conclusions 

We have examined the mid-IR and X-ray properties of 35 X-ray 
selected sources in the GOODS-N area which are optically faint 
(Rab > 26.5) and therefore missed in ground-based optical ob- 
servations. Instead of relying on previous work on the matching 
between X-ray and optical counterparts, our sample has been 
compiled anew. We find secure 3.6/mi counterparts for the X- 
ray sources using a likelihood ratio technique and then in turn 
we search for their possible optical ground-based Rab counter- 
parts. In case where there are none down to Rab = 26.5, we 
search for HST ACS optical counterparts. 18 sources have HST 
counterparts while the remaining have no optical counterparts. 
Our findings can be summarized as follows: 

1. Our sources populate moderate to high redshifts, being at 
markedly different redshifts from the already spectroscopi- 
cally identified population which peaks at z ~ 0.7. In par- 
ticular, the redshifts of the AGN with HST detections have 
moderate values with a median redshift of 2.6. The redshifts 
of the sources with IRAC detections only are definitely more 
uncertain; the objects with blue [5.8]-[8.0] colours are proba- 
bly located at redshifts comparable with the HST population, 
Z ~ 2.5, while the remaining sources could lie at z > 3.2. 
A couple of V a nd i dropouts exist in our sample (previ- 
ously reported by Lehme ret al.ll2005t) . we however propose 
that they are moderate redshift (z = 2 — 3) dust-extinguished 
AGN, rather than lying at very high redshift (z > 5). 

2. The sources with no optical counterparts in deep ground- 
based optical surveys constitute a large fraction (> 50%) of 
the total source population at high redshift (z > 2). This has 
important implications for the calculation and modelling of 
the luminosity function at high redshift, which in the case of 
our highly complete sample (97% of the X-ray sources have 
spectroscopic or photometric redshifts) is better represented 
by a modified PLE model. 

3. Our sources present very red colours. In particular, all 35 
sources with available Ks magnitudes would be character- 
ized as EROs on the basis of their Vgo6 - #s colour. 

4. The mid-IR colours of the sources with HST counterparts lie 
outside the AGN wedge, having blue [5.8]-[8.0] colours, in 
a region occupied by z ~ 2 galaxies according to the galaxy 
templates. The majority Spitzer IRAC mid-IR colours of the 
remaining sources with no optical counterparts lie within the 
AGN 'wedge', suggesting either QSO templates or galaxy 
templates at high redshift (z > 3). 

5. We find 4 high redshift candidates based on their non detec- 
tion with the HST, red [5.8] - [8.0] colour, and non detection 
with MIPS at 24 pm. However, the low IRAC 3.6 pm of 3 of 
them do not definately rule out dust obscuration in optical 
wavelengths. 

6. 17 out of 35 sources are detected in 24 /mi and can be clas- 
sified as optically faint mid-IR bright galaxies. This class of 
objects is widely believed to consist of reddened sources at 
moderate (~ 2) redshifts. 

7. The mean X-ray spectrum of our sources is very hard with 
T w 0.9, much harder than the spectrum of all sources 
in the CDFs suggesting that we are viewing heavily ob- 
scured sources. The X-ray spectroscopy on the individual 
sources suggests that three sources are candidate Compton- 
thick AGN. 

Obviously the current deepest X-ray observations are not 
at par with the present day optical spectroscopic capabilities. 
Spitzer has detected the faintest X-ray sources and thus provided 



aid in the determination of their properties and photometric red- 
shifts. 
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Table 1. List of CDFN X-ray sources which are correlated with a Spitzer 3.6yum detection and lack an optical counterpart in Capak et ^Tl J2004b . 
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62.262653 








>27.800 


22.156 


22.282 


9.228 


13.102 


19.188 


16.181 


<56.234 


>0.89 


4.90 m 


374 


1.69e-15 


189.307814 


62.307393 








>27.800 




>22.069 


3.490 


6.133 


9.344 


15.957 


55.452 


>1.47 


4.27 | 
2.44 2 


417 


7.35e-16 


189.356457 


62.285502 


27.678 


27.280 


26.292 


25.779 


22.512 


20.886 


12.484 


15.920 


20.284 


19.973 


172.580 


0.90 


434 


8.21e-16 


189.392075 


62.273511 


28.962 


27.275 


26.012 


25.042 


21.725 


20.434 


12.114 


12.983 


9.846 


6.153 


<56.234 


0.94 


3.39 8 


445 


4.43e-15 


189.404382 


62.309580 


>29.544 


28.945 


28.249 


27.339 


20.903 


21.146 


11.426 


14.165 


17.488 


15.500 


68.500 


2.34 


2.52 | 


470 


1.95e-15 


189.459284 


62.233180 


27.843 


27.448 


25.955 


24.867 


21.549 


20.254 


16.937 


19.801 


20.448 


21.565 


<56.234 


1.39 


1.67 I. 



